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ABSTRACT 


The objective of this study is to investigate the utility of acoustic tomography 
for performance assessment of a generic low frequency active sonar system. The 
performance of the sonar is simulated using tomography-derived sound speed data 
versus a range independent ocean model. The ocean environment used in the 
simulation is 159 tomographic snapshots of the Barents Sea Polar Front, taken 
every 5 minutes in August 1992. The modeled sonar system consists of a 1000 Hz 
source with a source level of 205.5 dB and a towed horizontal array of 
hydrophones. The system is derived from unclassified parameters of ATAS (Active 
Towed Array Sonar), built by Thomson Sintra ASM and British Aerospace SEMA, 
and the experimental ALF sonar, designed by FEL-TNO (the Netherlands) and built 
by Thomson Sintra ASM. The tomographic images over a range of 26 km provide 
a realistic ocean in which system performance is assessed. This study used a 
broadband, coupled normal mode, propagation model and assumed a noise-limited 
condition. The probability of detection calculated as a function of time for 13 hours 
is compared with that estimated using a range- and time-independent assumption. 


The utility of coastal acoustic tomography for tactical applications is discussed. 
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I. INTRODUCTION 


Since its introduction by Munk and Wunsch (1979), ocean 
acoustic tomography has evolved from the experimental phase 
into a phase in which different applications are feasible 
(Miller and Franken, 1994). Besides being of significant 
interest to oceanographic institutes, acoustic tomography can 
have applications for the world’s navies. They need range 
dependent propagation prediction models to assess the sonar 
performance in ocean areas with complicated Sound Speed 
Profiles (SSP). The increased importance of anti-submarine 
warfare in coastal ocean regions of the world requires more 
accurate and timely sound speed maps. Acoustic tomography can 
be utilized to provide space and time dependent maps of the 
ocean sound speed field. 

In order to quantify the above need for range and time 
dependent SSP’s, sonar performance assessment using acoustic 
tomography derived SSP’s is compared to the assessment using 
a conventional, range independent method, e.g. the use of 
expendable bathythermograph (XBT) data. The objective of this 
study is therefore to compare the sonar performance assessment 
of a low frequency active sonar system using tomography- 
derived sound speed data versus a range independent ocean 
model in order to show the utility for acoustic tomography as 
a tool for long range sonar performance evaluation. 

The tomographic and XBT ocean models are compared. The 
ecean environment used in the simulation is 159 tomographic 
snapshots of the Barents Sea Polar Front (every 5 minutes) 
taken in August 1992. The fixed location of the experiment is 
the Barents Sea Polar Front about 50 nm east of Bear Island 
(Figure 1). In the study, a 1000 Hz active sonar system with 
an acoustic power of 1 kW is assumed. The sonar performance 
assessment is quantified by the probability of detection on a 


LONGITUDE. 

Figure 1. Location of the Barents Sea 
Polar Front Experiment by the NPS, WHOI 
and SAIC in August 1992 (Miller, J.H. et 
al., 1994). 


submarine with a target strength of 10 dB. 

First the environmental setting of this study is 
described in Chapter II, "The Imaged Ocean." In Chapter III 
the hypothetical sonar system used for the comparison of the 
models is described. All parameters of the sonar are derived 
from unclassified data and, though the numbers are fictitious, 
they represent realistic figures for a modern low frequency 
active sonar system. Chapter IV deals with the transmission 
loss calculations by a coupled mode approach. The transmission 
loss for both, the range dependent and range independent, 
models is calculated using the same code. Finally the 
resulting probability of detection for noise-limited 
conditions is presented in Chapter v, followed by some 
concluding remarks in Chapter VI. 


II. THE IMAGED OCEAN 
A. WHY A RANGE DEPENDENT MODEL? 


In order to calculate the sonar performance assessment, 
we have to choose an appropriate acoustic model that is valid 
for the area of operation and requires a reasonable amount of 
computer calculation time. Since it is impossible in a real 
ocean to calculate exactly the real parameters, any model is 
an estimate. In general, range independent models are simpler, 
require less computing time, but are less accurate in range 
varying conditions. 

When we deal with range and time varying oceanographic 
conditions, we prefer a model that includes the variability of 
the conditions as a function of range and for which the 
calculations can be performed within the characteristic time 
slot of the changing conditions. In this case we need a range 
dependent model, that requires regular inputs of the changing 
conditions and fast calculations. 

Standard routine on board naval ships is to measure the 
acoustic conditions locally with an XBT or sound velocity 
meter (XSV). Usually, this is done in the operation area in 
six hour intervals or when the acoustic conditions may have 
changed substantially. We can use the resulting SSP as input 
for a range independent model. This method does not take into 
account the range variability of the ocean when we deploy long 
range sonar systems. For short range sonar though, this method 
can be adequate, since the range variability of the ocean is 
often negligible for these short ranges. We also can use the 
SSP’s of several ships at different locations as the input for 
a range dependent model. But still this method does only use 
spot measurements and does not take the time variability of 
the ocean into account. 

Sonar predictions for a long range sonar system in an 
ocean with varying acoustic conditions, e.g. near an 


oceanographic front, require a range dependent model with 
inputs of sufficient range and time resolution. These inputs 
can be provided by ocean acoustic tomography. 


B. OCEAN ACOUSTIC TOMOGRAPHY 


Ocean acoustic tomography is an efficient method to 
measure current or sound speed fields over a large ocean 
volume, described by Miller and Franken (1994), Spindel and 
Worcester (1990), and many others. It is based on travel time 
perturbations between sources and receivers of acoustic 
signals traveling along different paths. 

Acoustic tomography uses few transmitters and receivers 
to form a network of crossing raypaths in a horizontal area 
(Figure 2). By transmitting acoustic signals from one 
transmitter to any receiver over multiple paths in the 
vertical plane, one can deduce the properties of the ocean's 
interior in this plane on the basis of how the ocean altered 


the signals. Combining the information of all vertical slices 


Figure 2. The tomographic monitoring 
system: transducers (S) and receivers 
(R). Addition of one receiver (R’) 


increases the number of acoustic paths 
multiple times. 


of the network enables us to derive a three dimensional image, 
much like the medical CT scan. 

Adding more transducers or receivers increases the number 
of slices at a rate much greater than achieved with 
conventional instruments that do not collect data from 
crossing paths. By using transceivers, even more paths can be 
used and two way measurements can be done, enabling 
determination of current fields. Moreover these measurements 
are performed at the speed of sound, about 3000 knots in 
water, much faster than any research vessel can sail. It 
enables us to obtain a nearly continuous update of the ocean 
structure. 

Acoustic tomography can provide continuous, nearly real 
time information to predict sonar performance and counter 
detection ranges in a fixed area. Sometimes naval forces are 
required to operate in one area for a prolonged period of 
time. In these circumstances task group commanders need 
continuous information of the changing environment in this 
unchanging area for their tactical decisions. Detailed 
knowledge of the acoustic propagation in this area will be a 
key to success in the prosecution of submarines and in the 
hunting of mines. Future tomography instruments may be dropped 
by maritime patrol aircraft to monitor an area of interest. 
Also, with the upcoming interest in low frequency active and 
bistatic sonars (Delmee et al., 1994), one can think of using 
their signals for acoustic tomography, thus providing an 
onboard capability to monitor the operational area of a task 
group. Requirements for moving ship tomography are within 
current technical possibilities. They include very accurate 
position- and timekeeping, already available through GPS. For 
coastal acoustic tomography though, we need a vertical 
hydrophone array as receiver since we have to spatially 
beamform as the temporal processing of the data is not 
sufficient to identify the different signal arrivals (Miller 


et al., 1994). A vertical array is less easily usable froma 
surface vessel. 


BARENTS SEA POLAR FRONT EXPERIMENT DATA 


The Barents Sea Polar Front Experiment (BSPFEX), 
conducted by a team of the Naval Postgraduate School (NPS), 
Woods Hole Oceanographic Institution (WHOI) and Science 
Applications International Corporation (SAIC) in August 1992, 
is described by McLaughlin (1993). It provides this study with 
159 two-dimensional tomography snapshots (Figure 3), with a 
five minutes interval, of a coastal, range varying 
environment. When these snapshots are used as sequential 
frames for a movie, we observe a wave motion in the SSP of the 
frontal system, adding the temporal variation of the SSP to 
the study with a resolution of 5 minutes. The spatial 
resolution of the sound speed data is 1000 meters in range and 
2 meters in depth. 

From these tomographic snapshots, the XBT data are 
extracted by using the SSP of one location, extrapolated over 
the whole area (Figure 4). We have chosen the location of the 
source for this study at position "23 km" in the tomographic 
slice to be also the location where the XBT is taken. The XBTs 
are taken at the beginning of the data sequence and at the 
synoptical hours 00.00 and 06.00 GMT. This establishes the 
temporal differences between the tomography snapshots and the 
XBT based calculations. By comparing Figures 3 and 4, we can 
expect a substantial difference for the acoustical conditions 
downslope of this location, towards the frontal system, while 
upslope the differences between the tomography and XBT data 
are more subtle. 


NPS/WHOI Barents Sea Tomographic Image v2.0 


Depth (m) 
3 


1992 228: 20:02 UTC 


1470 


300 
0) 


Range (km) 
Figure 3. Tomographic snapshot. 


 NPSWHOI Barents Sea XBT image v2.0, position 23 


£20: 02 UTC 


1450 1460 1470 
Sound Speed (m/s) 
10 15 20 25 30 35 
Range (km) 
Figure 4, XBT image of the two-dimensional SSP 
derived from the tomographic snapshot. 


At this location there are three different water masses 
to distinguish: 


1. The warmer, shallow, mixed surface layer, extending to 
a depth of 20 meters. 


2. The cold, Arctic water, sound channel between 20 and 
120 meters depth. 


3. The North Atlantic water layer at the bottom of the 
coastal ocean. 

We studied three cases with the sonar lowered in one of above 
mentioned water masses. The used sonar depths are 10, 70 and 
180 meter as indicated in Figures 3 and 4 by the "@." For this 
last sonar depth we keep a bottom clearance of 45 meters. 

The bottom of the Barents Sea at this location consists 
of rock, covered by a 5 meter thick sediment. For the 
calculations we used the sound speeds and densities as 
displayed in Table 1. 

The studied sonar range extended 13 kilometers up- and 
down slope from the sonar position as indicated in Figures 3 
and 4 by the dashed line. 


sound speed (m/s) | density (kg/m’) 
water i.a.w Figure 3 1026 
sediment 1659 1830 
bottom 4213 2600 


Table 1. Sound speed and density data for the 
BSPFEX (sediment and bottom data obtained from US 
Naval Oceanographic Office: NAVOCEANO) . 


D. THE ENVIRONMENTAL ASSUMPTIONS 


In order to calculate the transmission loss as a function 
of range and depth, we applied absorption to the pulse 
propagation. The absorption figures are derived from Clay and 
Medwin (1977). The modal attenuation coefficient is calculated 
for the frequency of 1000 Hz by the formula: 


a, 
8, = [pep Bz) a2) ae ; q@ 
where Z,(z) is calculated by the normal mode code, @ = 0.3 
dB/m in the sediment and @ = 8-10°°dB/m in seawater. The depth 
of the waveguide, h in Equation (1), is for this study 300 
meters. 

Finally the modal attenuation coefficient is corrected 
for bottom and surface interactions. The coherent reflection 
coefficient for a surface having a Gaussian PDF is given by 
the following Equation (Clay and Medwin, 1977): 


R, = Re 2eosin'e, (2) 


where R = 1 for surface interactions and R = 0.759 for bottom 
interactions. For the Barents Sea o = 0.25 meter for the 
bottom. For the surface we chose 6 = 0.5 meter, corresponding 
with sea state 3. The bounce distance for each mode is 


calculated by (Chiu, 1994): 


with which the number of bounces for each mode over the range 
of 13 km is calculated. 

For a source depth of 70 meters, these corrections 
contribute to the attenuation coefficient for mode 18 and 
higher modes. As shown in Figure 5, the attenuation for modes 
higher than 75 is larger than 0.5 Np/km. This is one reason to 


limit the number of modes for the transmission loss 
calculations to 75. 

For the probability of detection calculations we use a 
noise specrum level, NSL = 60 GB re 1pPa/Hz*, as derived from 
the Wenz-curve (Kinsler et al., 1982) for sea state 3. 


10° ‘modal attenuation coefficient for 1000 Hz 
/ 
3rincluding bottomvsurtace losses: 
‘sigma(surtace) 5m (SS 3) 
és sigma(bottom ) = 0.25 m 
without bottom/surtace losses: 
= 
Hy 2 
2 
15 
1 
os} ~~ -------------=5 
° 
oO 50 75 100 180) 
modo# 
Figure 5. Modal attenuation coefficients 
for 1000 Hz. 
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III. THE HYPOTHETICAL SONAR SYSTEM 
A. LOW FREQUENCY ACTIVE SONAR 


Today’s extremely quiet nuclear submarines as well as 
widespread conventional submarines are requiring many navies 
in the world to have research programs for low frequency 
active sonar systems. Threat scenarios have changed to area 
operations in coastal regions with a subsurface threat of very 
modern to old conventional submarines and of mines (O'Keefe, 
1992). Current sonar systems do not provide sufficient 
detection capabilities in these scenarios. The submarines are 
too quiet for passive detection and their range advantage over 
ships fitted with hull mounted sonars is increasing. 

A low frequency active sonar, with frequencies below 1000 
Hz, can provide long range detection capabilities against 
quiet submarines, primarily in deep water. It can be designed 
monostatic, with source and receiver in one unit, or 
multistatic, with separate units for source and receivers. In 
confined and shallow waters a better option is an activated 
towed array sonar (ATAS), operating at slightly higher 
frequencies. Such a system by Thomson Sintra ASM and British 
Aerospace SEMA is currently on the market, operating at a 
frequency around 3 kHz (Jane’s, 1994). The use of a low 
frequency active sonar in shallow water is under investigation 
(Delmee et al., 1994). 


B. SONAR SYSTEM DESIGN FOR THIS STUDY 


The bases of the sonar system design for this study are: 


1. The experimental ALF sonar, designed by FEL-TNO in the 
Netherlands and built by Thomson Sintra ASM. The 
system is described by De Vlieger et al. (1994). 
Currently the Dutch Navy is testing the ALF in sea 
trials, planned from 1994 until 1997. 
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2. ATAS, designed and built by Thomson Sintra ASM and 
British Aerospace SEMA. The parameters are extracted 
from an undated sales brochure and Jane’s (1994). 


3. Suggestions of Directorate Material Royal Netherlands 
Navy, Department of Weapon and Communication Systems, 
division OSAO. 


The system consists of a 1000 Hz towed source and, attached to 
the source body, a single receiver array. The design is 
similar to ATAS, using a hypothetical frequency. The 
parameters of the source are derived from the ALF sonar and 
shown in Table 2. 


| Frequency 1000 Hz 

| Source acoustic power 1 kw 

| horizontal beampattern omni-directional | 
vertical beamwidth 30° 
DI 4.5 dB 
sL 205.5 dB 

| array dimensions 0.84x2.50 m 

| # of transducers 2x6 
pulse FM up-chirp 
pulse length 10 sec 
pulse swept bandwidth 100 Hz 

Receiver horizontal beamwidth 2° 

AG 15 dB 
array length 21.7 m 
# of hydrophones 32 


Table 2. The hypothetical sonar system parameters. 
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1. The Source 


The source consists of two vertical line arrays of six 
tonpilz transducers with a length of 2.5 meters as shown in 
Figure 6. It transmits a CW pulse at f£ = 1000 Hz with an 
acoustic power of 1 kW. For such a planar array, assuming no 
complex weighing, the elements to be omni-directional and for 
a speed of sound c = 1450 m/s, the far-field directivity 
function is: 


3 


2 
Dif. f,£) => 
= 


e2R (1m 0.5) Laer (2 0.51 fae) (4) 
4 


where: 


£, = sin(8) cos(w) = , (5) 


£,= cos (0) t (6) 
é 


and the spacing between the elements dx=dz=42 cm. @ is the 
angle describing elevation and w represents the azimuthal 
angle. From the plot of the directivity function as shown in 
Figure 7, we can estimate the directivity index (Kinsler et 
al., 1982): 


T= L0Log.(F-) + (7) 


where 6’ is the -6dB half beam width in radians. The 
directivity index for the sonar, DI = 4.5 dB. The resulting 
source level, SL = 205.5 dB. 

From the same figure, we measure the -3dB beamwidth to be 
30°. To minimize the calculation time of the transmission loss 
calculations, we choose the half beamwidth to determine the 
maximum number of used normal modes. The angle of propagation 
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tow cable 


pressure 
box 


tow points 

‘ towed arrays 
Tonpilz 

transducer 


Figure 6. The source as copied from the experimental 
ALF sonar (De Vlieger et al., 1994). 
(© Thomson Sintra ASM, reproduced by permission) 


ALF source vertical beampattern for 1000 Hz 


20) 


g 
= 
40) 
-s0} broadside (__) endfire(..) 
beamwidth = 28 dgr beamwidth = 20 dgr 
Di=4608 DI=4.408 | 
i a re a 
theta (elevation) 
Figure 7. The directivity function for the 
hypothetical sonar source for a frequency of 
1000 Hz. 
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of mode# 75, 6, = 14.9°, thus providing a second reason to 
limit the calculations to 75 modes. 

Use of low frequencies requires relative long pulses to 
meet a certain detection threshold. In order to have a 
reasonable range and doppler resolution, most low frequency 
active sonar systems will use frequency (FM) or phase (PM) 
modulated pulses. The hypothetical sonar system uses an up 
chirping, linear FM pulse with a 100 Hz swept bandwidth of 
length 10 seconds. But for the transmission loss calculations 
of this study, we modelled the source as a point source and 
limited the pulse to be a CW pulse at one frequency, again to 
minimize the computer calculation time. 


2. The Receiver Array 


The single receiver array with a length of 21.70 meters 
consists of 32 hydrophones with an elememt spacing of 
approximately half a wave length. It is towed behind the 
source as indicated in Figure 6. As shown in Figure 8, the 
broadside horizontal beamwidth is 2 degrees and the array 
gain, AG = 15 4B. 


15 


ALF receiver array broadside beampattemn 


; 
frequency = 1000 Hz 
temaeesae | 

ng array length = 21.7m_ 

20) | 


60 80 100 120 140 160 
psi (azimuth) 


Figure 8. The directivity function for the 
hypothetical sonar receiver. 
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Iv. THE TRANSMISSION LOSS CALCULATIONS 


In order to assess the sonar performance in the ocean, we 
have to model the sound propagation through an inhomogeneous 
medium varying in four dimensions: three spatial variables and 
one time variable. In the imaged ocean of this study, the 
tomographic image, the soundspeed field varies in range, depth 
and time. We take care of the temporal fluctuations by 
calculating 159 sequential snapshots of the ocean with a 5 
minutes spacing. Still the remaining two varying dimensions, 
range and depth, require a range dependent model for the sound 
propagation calculations. In case of shallow water and low 
frequency, a normal mode approach is suitable. Although for 
normal mode theory 1000 Hz is not considered to be a low 
frequency, the coupled normal mode model is particularly 
useful for active sonar applications (Chiu et al., 1994) as 
will be explained next. 


A. COUPLED MODE APPROACH 


The imaged ocean can be seen as a waveguide with 
significant environmental range variation; the SSP’s are 
highly range dependent, especially across the front, and the 
bathymetry is varying. In such a case a travelling mode in the 
waveguide exchanges energy with other modes as explained by 
Chiu (1994). When it travels downslope over a gentle slope, 
the mode will exchange energy with the mode of one order 
lower. The modes are coupled. This requires a range dependent 
propagation model that calculates the propagating modes as a 
function of range. 

The newly developed 3-D method we use in this study is a 
three-dimensional, broadband, coupled normal mode sound 
propagation model, described by Chiu et al. (1994) and 
provided to me by EOS Research Associates of Monterey, 
California. It assumes the ocean to be a linear, time- 
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invariant system that filters the input source signal and 
produces the output three-dimensional sound field. The model 
determines the ocean transfer function at various acoustic 
frequencies. To cut down in computation time we use a CW 
source at one frequency, 1000 Hz. 


For a given frequency, the basic formulation 
involves decomposing the acoustic pressure into 
slowly varying complex envelopes that modulate 
(mode by mode) analytic, rapidly-varying, 
adiabatic-mode solutions.... Since the formulation 
allows for sound speed and density to vary three- 
dimensionally and bathymetry to vary two- 
dimensionally, it is particular useful for...active 
sonar applications. (Chiu et al., 1994) 


Modal attenuation coefficients are implemented to account for 
sound absorption in seawater and sediments, and boundary 
scattering (see Chapter II). Backscattering is neglected. 

When only a small number of modes are required for the 
calculation, i.e. in shallow water and for low frequencies, 
the method is very effective. Since the frequency 1000 Hz 
stretched the model to its computational limits, we have 
reduced the number of used modes to 75 as justified in Chapter 
II and III. With the limitations imposed on the computations 
for this study, it takes on average 1°08’ per image to 
calculate the transmission loss for each frame, using a HP 
Apollo, series 735, workstation. While this is much longer 
than the 5 minutes interval between the subsequent tomographic 
data measurements, developments in computing hardware, 
software, and refinements in the propagation code could make 
this technique feasible in real time. 


B. SET-UP OF THE TRANSMISSION LOSS COMPUTATIONS 


We developed a code to calculate the transmission loss in 
a waveguide representing the imaged ocean of Figure 3. It 
sequentially calculates transmission loss for each snapshot, 
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including the XBT images at the beginning of the sequence and 
at the standard synoptical hours. We use the same coupled 
normal mode model to calculate the transmission loss for the 
tomographic and XBT snapshots. This implies that the 
calculations for the XBT snapshots are not range independent 
since the model takes the varying bathymetry into account. But 
this is the only fair methodology to compare both the 
tomographic and XBT images. Moreover the bathymetry data is 
normally available on board naval ships. 

The computations are performed twice in a waveguide of 
depth 300 meters and length 13 kilometers: starting from the 
source, once downslope and once upslope as indicated in Figure 
3. The input of the code consists of the source position in 
the imaged ocean and the sound speed and density fields in the 
water, sediment and bottom with a range increment of 1000 
meters and a depth increment of 2 meters. For the computations 
as represented in this thesis, all other parameters to the 
coupled normal mode model are kept constant. 

The broadband coupled normal mode model numerically 
integrates the differential equations governing the complex 
envelopes that describe the acoustic pressure, in the radial 
direction. To allow for economical computations, the range 
step size is automatically adjusted to the horizontal 
environmental variations. We choose the depth increment to be 
0.75 meters, approximately half of the wavelength. This was 
the maximum possible depth step for calculations up to mode 
75. Plots of the local normal modes showed aliasing effects 
for mode 100 and higher modes as can be seen in Figure 9f. 
Figures 9a-f show another interesting feature of the local 
normal modes close to the source. Low modes are more affected 
by the SSP while high modes are more affected by the 
bathymetry. 

For each frequency, the code calculates the farfield 
sound pressure as a function of range, r, and depth, z: 
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Figure 9f. Local normal 
mode for mode 120. 


pir.z) = PP or) Z,(z.0) SO (8) 
yom 

where P, is the rms pressure at 1 meter from the source, N is 
the number of propagating modes, U,’s are the slowly varying 
complex envelopes of the horizontal structure of the modes, 
Z,!S are the local normal modes and k,’s are the modal 
horizontal wave numbers. These variables are all given by the 
broadband coupled normal mode model. Hence the transmission 
loss, TL, can be calculated, using: 


TL(r,z) = ~10log,,{pir,z) p*(r,z)) (9) 


where * denotes the complex conjugate. The acoustic pressure 
is a complex quantity since we include attenuation. 


. RESULTS 


The output, the transmission loss, is plotted with a 
range increment of 500 meter and a depth increment of 2 meter. 
Two-hourly plots are shown in Appendices A, B and C for the 
respective source depth of 70, 180 and 10 meters. Additionally 
a movie of the 5-minute plots was made (Franken, 1994). 

An example of the temporal transmission loss variations 
for a fixed position (range and depth) in the imaged ocean is 
shown in Figure 10. Besides the spatial difference in 
transmission loss between the tomographic and XBT calculations 
based on the choice of the location, Figure 10 shows clearly 
the temporal fluctuations of the tomography derived 
transmission loss. These fluctuations are wiped out for the 
XBT calculations since the XBT is only launched once every six 
hours. 

In order to get a real flavor of the temporal 
transmission loss variations, it is more interesting to watch 
the movie mentioned above. While the movie of the SSP displays 
clearly a wave motion in the front of cold and warm water, the 
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‘Temporal TL Variation for Selected Resolution Cell 


TL (aB) 


oO = 4 6 8 10 12 
time (hrs) 
Selected resolution cell: tomography = 
range: 8 km xeT 
depth: 130 m 
Figure 10. The transmission loss temporal 


variability at position: range 8 km downslope, 
depth 130 m. 


transmission loss movie shows a more erratic behavior, 
especially downslope from the source towards the front. The 
transmission loss at a certain range from the source is 
affected cumulatively by all variations in the sound speed 
profiles between that range and the source. Many sound speed 
variations between source and field point makes’ the 
transmission loss vary randomly. Comparing subsequent 
snapshots, a temporal variation is clearly visible. 
Comparing the XBT and tomography snapshots of the 
transmission loss, two obvious observations can be made. 


1. In the tomographic snapshots the effects of the front 
are clearly visible: the energy of the sound refracts 
downward towards the bottom where the sound speed is 
the lowest. A spatial difference between the two 
snapshots in downslope direction can be observed. 
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2. Although the sound speed fields of the XBT and 
tomographic images do differ in range and time upslope 
from the source, the general structure of both sound 
fields remains the same and thus the transmission loss 
snapshots are similar. A spatial difference between 
the two snapshots in upslope direction is not clear. 


These two observations are visualized in Figure 11. This 
Figure shows the mean difference in transmission loss between 
the tomographic and XBT calculations over the 159 available 
snapshots. The mean difference in transmission loss is 
calculated by the formula: 

2 


Pri 
Thaser = ieee ee , (10) 


‘7 


where P,,? and Py? are the acoustic intensities for the 
tomographic and XBT calculations respectively as a function of 
range and depth for snapshot "i." 

We thus conclude that if the sound speed field does not 
change significantly in time and range, the XBT based 
transmission loss calculations and subsequent sonar assessment 
are similar to the tomographic snapshots. But in case of 
drastic sound speed field changes, in time or range, acoustic 
tomography will contribute a lot to a more accurate, range 
dependent transmission loss calculation. This is not only 
important in the case of this study, where a front changes the 
sound speed field drastically, but also for long range sonars 
in open ocean, where the sound speed structure of the sea 
changes gradually, but over a long distance from the source to 
its maximum range. It is also very important when temporal 
changes take place, unnoticed by the 6-hourly XBT, e.g. in 
case of a afternoon effect. 
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Figure 11. Mean difference 
transmission loss between the tomogr: 
and XBT calculations over 159 snapshots. 


V. THE SONAR PERFORMANCE 


An operational measure of sonar performance is the 
probability of detection. The probability of detection is 
defined as the probability of detecting a signal in the 
presence of noise, also called the single ping probability of 
detection or the glimpse probability. Clearly this is not the 
probability of target detection, since that may be a function 
of the probabilities for several pulses. Target detection 
depends also on many outside factors like relative movement of 
the target, sonar operator performance, etc. For sonar system 
performance assessment though, the single ping probability of 
detection is most appropriate. For the detection of the signal 
in the presence of noise we do not take any human factor of 
the sonar operator into account. 

The single ping probability of detection can be derived 
using the concept of hypothesis testing as explained by Burdic 
(1984). Using the Neyman-Pearson criterion, we can express the 
probability of detection, P,, as a function of the probability 
of false alarm and the signal-to-noise ratio (SNR) (Ziomek, 
1985): 


P pam (11) 


The probability of false alarm used in this study, Py, = 10°. 
The SNR is derived from the noise-limited (monostatic) active 


sonar Equation: 
SNR(dB) = SL-2TL+TS-NL-DI, (12) 
where the source level, SL, is (Kinsler et al., 1982): 


SL{relpPa) = 10logII+DIr+171, (13) 


for the acoustic power, Il = 1 kW, and the source directivity 
index, DI = 4.5 dB, as calculated in Equation (7). The 
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transmission loss, TL, is obtained by the coupled normal mode 
model. The target strength is a function of many target 
parameters and the frequency. A target strength, TS = 10 dB, 
is a reasonable assumption for a 60 meters long submarine, 
pinged abeam at 1 kHz (De Vlieger et al., 1994). The noise 
level, NL (see Chapter II), is calculated for the frequency 
bandwidth of 100 Hz. The directivity index, DI, in Equation 
(11) is the array gain, AG = 10log,,(32 elem.) = 15 4B. 

Equation (11) assumes noise limited conditions. We have 
not looked into reverberation limited conditions. This does 
not limit the research as long as we compare the tomographic 
and XBT snapshot under the same conditions. 

Appendices A to C show the probability of detection for 
the source depths 70, 180 and 10 meters respectively. When the 
spatial and temporal variability of the ocean are large, in 
the downslope region from the source, there is a large 
difference in the probability of detection between the 
tomographic and XBT snapshots visible. Under these 
circumstances it is favorable to use tomography sound speed 
data for the sonar performance assessment. 
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VI. CONCLUSIONS 


When the ocean sonar conditions vary significantly in 
range and time, a more realistic sonar performance assessment 
is to be expected using tomography sound speed data instead of 
XBT data. This becomes more important when long range sonars, 
such as low frequency active sonar, are involved. The ocean 
conditions are liable to change over the long maximum range of 
such a sonar system. This is also very important when temporal 
changes of the ocean conditions take place. An example of such 
a case besides the frontal system we use in this study, is the 
afternoon effect that can effect the performance of some 
sonars on a large scale. 

The comparison between the sonar performance assessment 
using tomography data and XBT data is fair since we use the 
same model and assumptions for both. This study is limited 
though: we look only at the sonar performance in one scenario 
under noise limited conditions. The differentiation in source 
depth does not contribute much to the conclusions. For the 
source at 10 meter, the maximum sonar range is poor and thus 
the ocean conditions do not change significantly over this 
small range. 

There will be a larger difference between the XBT SSP 
derived from tomographic measurements and from a real XBT- 
reading than assumed in this study. The sound speed does not 
solely depend on the ocean temperature, thus introducing a 
systematic error in the XBT derived sound speed data. This can 
be avoided by using an XSV measurement. 

Comparing the transmission loss computer computation time 
(1508’) and the temporal resolution of the tomography 
snapshots (5 minutes), it is clear that the sonar performance 
assessment for this 1000 Hz sonar using this broadband, 
coupled normal mode model, can not be performed real time. 
Though for lower frequencies, requiring less normal modes for 
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APPENDIX A. TL AND POD ESTIMATES FOR SOURCE DEPTH 70 METER 
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APPENDIX B. TL AND POD ESTIMATES FOR SOURCE DEPTH 180 METER 
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APPENDIX C. TL AND POD ESTIMATES FOR SOURCE DEPTH 10 METER 
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